Time parameters and Lorentz transformations of relativistic stochastic processes 
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Rules for the transformation of time parameters in relativistic Langevin equations are derived and 
discussed. In particular, it is shown that, if a coordinate-time parameterized process approaches 
the relativistic Jiittner-Maxwell distribution, the associated proper-time parameterized process con- 
verges to a modified momentum distribution, differing by a factor proportional to the inverse energy. 
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Stochastic processes (SPes) present an ubiquitous tool 
for modelling complex phenomena in physics 11], 0, Q , bi- 
ology [4, 5] , or economics and finance [a, 0, H, U Stochas- 
tic concepts provide a promising alternative to determin- 
istic models whenever the underlying microscopic dy- 
namics of a relevant observable is not known exactly 
but plausible assumptions about the underlying statis- 
tics can be made. A specific area where the formu- 
lation of consistent microscopic interaction models be- 
comes difficult [13, EH, m 

concerns classical relativistic 
many-particle systems. Accordingly, SPes provide a use- 
ful phenomenological approach to describing, e.g., the 
interaction of a relativistic particle with a fluctuating en- 
vironment 0, BE Bill]- Applications of stochastic 
concepts to relativistic problems include thermalization 
processes in quark-gluon plasmas, as p roduced in rela- 
tivistic heavy ion colliders [3, [3, [3. 21 L or complex 



high-energy processes in astrophysics [22l. l23l. 124 l25|. 

While these applications illustrate the practical rele- 
vance of relativistic SPes, there still exist severe con- 
ceptual issues which need clarification from a theoretical 
point of view. Among these is the choice of the time- 
parameter that quantifies the evolution of a relativistic 
SP [2^|. This problem does not arise within a nonrel- 
ativistic framework, since the Newtonian physics postu- 
lates the existence of a universal time which is the same 
for any inertial observer; thus, it is natural to formulate 
nonrelativistic SPes by making reference to this universal 
time. By contrast, in special relativity 27|, \2M the notion 
of time becomes frame-dependent, and it is necessary to 
carefully distinguish between different time parameters 
when constructing relativistic SPes. For example, if the 
random motion of a relativistic particle is described in 
a ^-parameterized form, where t is the time coordinate 
of some fixed inertial system E, then one may wonder 
if/how this process can be re-expressed in terms of the 
particle's proper-time r, and vice versa. Another closely 
related question 17j concerns the problem of how a cer- 
tain SP appears to a moving observer, i.e.: How does a 
SP behave under a Lorentz transformation? 

The present paper aims at clarifying the above ques- 



tions for a broad class of relativistic SPesgoverned by 
relativistic Langevin equations 3 3, 15, 16, 3] ■ First, 
we will discuss a heuristic approach that suffices for most 
practical calculations and clarifies the basic ideas. Sub- 
sequently, these heuristic arguments will be substanti- 
ated with a mathematically rigorous foundation by ap- 
plying theorems for the time-change of (local) martingale 
processes [29| . The main results can be summarized as 
follows: If a relativistic Langevin-Ito process has been 
specified in the inertial frame E and is parameterized 
by the associated E-coordinate time t, then this process 
can be reparameterized by its proper-time r and the re- 
sulting process is again of the Langevin-Ito type. Fur- 
thermore, the process can be Lorentz transformed to a 
moving frame E', yielding a Langevin-Ito process that is 
parameterized by the E'-time t' . In other words, similar 
to the case of purely deterministic relativistic equations 
of motions, one can choose freely between different time 
parameterizations in order to characterize these relativis- 
tic SPes - but the noise part needs to be transformed 
differently than the deterministic part. 

Notation.- We adopt the metric convention (r] a p) — 
diag(— 1, 1, . . . , 1) and units such that the speed of light 
c = 1 . Contra- variant space-time and energy- momentum 
four- vectors are denoted by (x a ) = (x°,x l ) = (x°,x) — 
(t,x) and (p a ) — (p°,p l ) — (p°,p), respectively, with 
Greek indices a — 0,1,..., d and Latin indices i = 
1, . . . , d, where d is the number of space dimensions. Ein- 
stein's summation convention is applied throughout. 

Relativistic Langevin equations.- As a starting point, 
we consider the ^-parameterized random motion of a 
relativistic particle (rest mass M) in the inertial lab 
frame E. The lab frame is defined by the property 
that the thermalized background medium (heat bath) 
causing the stochastic motion of the particle is at rest 
in E (on average). We assume that the particle's 
trajectory (X(t),P(t)) = {X l (t), P z (t)) in E is gov- 
erned by a stochastic differential equation (SDE) of the 
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form [IE El BSE Gil 



dX a (t) = {P a /P°)dt, (la) 
dP l (t) = A 1 ,11 • f, d/,":/ :. (lb) 



Here, dX°(t) = dt and dX l (t) := X*(t + dt) - X l (t) 
denote the time and position increments, dP l {t) 
P { (t + dt) - P*(i) the momentum change. P°(t) 
(M 2 + P 2 ) 1 / 2 i s the relativistic energy, and V*(i) 
dX l /dt = P l I P° are the velocity components in E. In 
general, the functions A 1 and C l j may depend on the 
time, position and momentum coordinates of the parti- 
cle. The random driving process B(t) — (B J (t)) is taken 
to be a d-dimensional ^-parameterized standard Wiener 
process 29|, 3(1 31| . i.e., B(t) has continuous paths, for 
s > t the increments are normally distributed, 

-\u\ 2 /[2 (s-t)] 

?{B(s) - B(t) e[u,u + du}} = ^ _ /2 d d u, (2) 



and independent for non-overlapping time intervals 4l| . 

Upon naively dividing Eq. (|lb[) by dt, we see that A 1 
can be interpreted as a deterministic force component, 
while C l jdB 3 (t)/dt represents random 'noise'. However, 
for the Wiener process the derivatives dB 3 (t) / dt are not 
well-defined mathematically so the differential represen- 
tation JT]) is in fact short hand for a stochastic integral 
equation [29|, [3l| with C l dB 3 signifying an infinitesimal 
increment of the Ito integral [32L |33| . Like a determin- 
istic integral, stochastic integrals can be approximated 
by Rlemann-Stieltjes sums but the coefficient functions 
need to be evaluated at the left end point t of any time 
interval [t, t + dt] in the Ito discretizati on [4211. In con- 
trast to other discretization rules [J 0, l3lL I34L l35| , the 
Ito discretization implies that the mean value of the noise 
vanishes, i.e., (C l jdB 3 (<)) = with ( • ) indicating an av- 
erage over all realizations of the Wiener process B(t). In 
other words, Ito integrals with respect to B(t) are (local) 
martingales [2{|. Upon applying Ito's formula [Il[3l| to 
the mass-shell condition P°(t) = (M 2 + P 2 ) 1 / 2 , one can 
derive from Eq. 
tivistic energy: 



rb|) the following equation for the rela- 



dP°(t) =A°dt + C° r dB r (t), 



(3) 
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where A, := A\ D i3 := D lj = C l C l an d dr ■= & r . 

Equations CO) define a straightforward relativistic 
generalization El J3j EH of the classical Ornstein- 
Uhlcnbeck process [36| , representing a standard model of 
Brownian motion theory [43j . The structure of Eq. (fTa|) 
ensures that the velocity remains bounded, |V| < 1, 
even if the momentum P were to become infinitely large. 
When studying SDEs of the type ([I]), one is typically in- 
terested in the probability f(t, x, p) d d xd d p of finding the 
particle at time t in the infinitesimal phase space interval 



[x, x + dx] x [p,p+dp\. Given Eqs. {T]), t ne non-negative, 
normalized probability density f(t,x,p) is governed by 
the Fokker-Planck equation (FPE) 



d ^ p l d \ _ d 
dt p° dx % J dp 1 



-A*f 



lJL(D ik f) 



(4) 



= where f i s a Lorentz scalar [33| and p° = (M 2 + 
= p 2 ) 1 / 2 !4J]. Deterministic initial data X(0) = xq and 
P(0) = p Q translates into the localized initial condition 
f(0,x,p) = S(x — xq) S(p — p ). Physical constraints on 
the coefficients A l {t, x,p) and C l r (t, x,p) may arise from 
symmetries and/or thermostatistical considerations. For 
example, neglecting additional external force fields and 
considering a heat bath that is stationary, isotropic and 
position independent in E, one is led to the ansatz 



A 1 = -a(p°) p l , C l .j = pD^ )] 1 / 2 S l 



(5a) 



where the friction and noise coefficients a and D depend 
on the energy p° only. Moreover, if the stationary mo- 
mentum distribution is expected to be a thermal Jiittncr 



function [38l. |39|. i.e., if foe := lim^oo / c>c exp(— (3p°) in 
S, then a and D must satisfy the fluctuation-dissipation 
condition M 



= a(p°) p° + dD{p°)/dp° - f3D{p°) 



(5b) 



In this case, one still has the freedom to adapt one of the 
two functions a or D. 

In the remainder, we shall discuss how the process ([T]) 
can be reparameterized in terms of its proper-time r , and 
how it transforms under the proper Lorentz group [28| . 

Proper-time parameterization.- The stochastic 
proper-time differential dr(i) = (1 — V 2 ) 1 ^ 2 dt may be 
expressed as 



dr(t) = (M/P°) dt. 



(6a) 



The inverse of the function r is denoted by X°(t) = t(r) 
and represents the time coordinate of the particle in the 
inertial frame E, parameterized by the proper time r. 
Our goal is to find SDEs for the reparameterized pro- 
cesses X a {r) := X a (t(r)) and P Q (r) = P a {t{r)) in E. 
The heuristic derivation is based on the relation 



dB j {t) 



p0\ 1/2 
"' =' M 



p0\ 1/2 



d& (r), (6b) 



where B 3 (r) is a standard Wiener process with time- 
parameter t. The rigorous justification of Eq. (|6b[) is 
given below. Inserting Eqs. ((6j in Eqs. |lj one finds 

dX a (r) = (P a /M)dT, (7a) 
dP*(r) = A 1 dT + C^d&ir), (7b) 

where A 1 := (P°/M) A l (X°, X, P) and C l j := 
(pv/M) 1 / 2 C l AX ,X,P). The FPE for the associated 
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Note that di'(i) = dY'°(t) = A° At dX' i (t), and, hence, 



FIG. 1: 'Stationary' probability density function (PDF) of 
the absolute momentum \P\ measured at time t = 15 (x) 
and t = 15 (o) from 10000 sample trajectories of the one- 
dimcnsional (d = 1) relativistic Ornstein-Uhlenbeck pro- 
cess [l3 |. corresponding to coefficients D(p°) = const and 
a(p°) = /3D/p° in Eqs. (fjQ) and ©. Simulation parameters: 
dt = 0.001, M = c = /3 = D = 1. 



probability density f(T,x°,x,p) reads 



d p a d 



ch 



M dx a 



f = 



d_ 

dp 1 



-A l f 



1 d 
28^ 



D ik f) 



(8) 



where now D lk := J^r^r^r- We note that 
f(r,x a ,x,p) dx°d d xd d p gives probability of finding the 
particle at proper-time r in the interval [t, t + dt] x [x, x + 
dx] x [p,p + dp] in the inertial frame E. 

Remarkably, if the coefficient functions satisfy the con- 
straints ([5]) - so that the stationary solution f^, of Eq. ((4]) 
is a Jiittner function 4>j(p) = exp(— (3p°) - then 
the stationary solution f^ of the corresponding proper- 
time FPE ([8]) is given by a modified Jiittner function 
4>mj(p) — exp(— (3p°)/p°. The latter can be derived 
from a relative entropy principle, using a Lorentz invari- 
ant reference measure in momentum space (4fj| . Physi- 
cally, the difference between and foo is due to the fact 
that measurements at t = const and r = const are non- 
equivalent even if r, t — ► oo. This can also be confirmed 
by direct numerical simulation of Eqs. ([I]), see Fig. [TJ 

Having discussed the proper-time reparameterization, 
we next show that a similar reasoning can be applied to 
transform the SDEs |T]) to a moving frame E' [l7| . 

Lorentz transformations.- Neglecting time-reversals, 
we consider a proper Lorentz transformation [28} from 
the lab frame E to E', mediated by a constant matrix 
A v ^ with A°o > 0, that leaves the metric tensor n a /3 
invariant. We proceed in two steps: First we define 

Y'"{t) := A\X»{t) , G lv {t) := A" M P"(t). 

Then we replace t by the coordinate time t' of E' to obtain 
processes X' a {t') = Y' a {t(t')) and P' a {t') = G' a (t{t')). 



dt'it) = ^ dt =^dt= P '° {t ' {t)) 



dt, (9) 



where A 1 is the inverse Lorentz transformation. Thus, 
a similar heuristics as in Eq. (|6b|) gives 



dB j (t) 



p o \i/2 
P*)^ 



p/0 



1/2 



where B'i{t') is a Wiener process with time t'. Further- 
more, defining primed coefficient functions in E' by 

A\ A" ((A- 1 ) ^, (A-%a/», (A- 1 ) 1 ^) , 
C')(x'°,x',p') := [(A- 1 ) ^"/^] 1 ^ x 
A^ ((A- 1 ) ^'^, (A^)>^, (A- 1 )*^") , 

the particle's trajectory (X'(t'),P'(t')) in E' is again 
governed by a SDE of the standard form 

dX la {t') = (P' a /P'°)dt', (11a) 
dP'\t') = A H dt' + C'\ J dB' i {t'). (lib) 

Rigorous justification.- We will now rigorously de- 
rive the transformations of SDEs under time changes and 
thereby show that the heuristic transformations leading 
to Eqs. and (fTTj) are justified; i.e., we are interested 
in a time-change 1 1— > t of a generic SDE 

dY(t) =Edt + Fj dB j (t), (12a) 

where E and Fj will typically be smooth functions of 
the state-variables (Y, . . .) HU, and B(t) = (B j (t)) is a 
d-dimensional standard Wiener process [46| . We con- 
sider a time-change t h- > i specified in the form [cf. 
Eqs. KBaJ) and ©] 



dt = H dt, 



t(0) = 0, 



(12b) 



with H being a strictly positive smooth function [47[ 
of (Y, . . .). The inverse of i(t) is denoted by t(t). We 
would like to show that Eq. (|12ap can be rewritten as 

d Y (i) = Edi + Fj dB j (t) , ( 1 2c) 

where Y{i) := Y(t(t)), E{i) := E(t(t))/ H(t(i)), 



(12d) 



F'(t) :=Fi(t(t))/^H(t(t)), and 

dB J '(i) :=VHdB j {t) 



is a (i-dimensional Wiener process with respect to the 
new time parameter i [48|. 

First, we need to prove that Eq. (|12d|) or, equiva- 

lently, := /*® dP J (s) does indeed define 



4 



a Wiener process. To this end, we note that for fixed 
j G {1, . . . , d} the process U (i) := J Q yjH(s) dB j (s) is a 
continuous local martingale, whose quadratic variation 



[L j ,L j ]{t) 



2"-l 

Um V <! V 

n — >og — * 

k=Q 



(k + l)t 



kt 



is given by [U,U](t) = f*H(s)ds @. For the 
quadratic variation of &{t) = LP(t(t)) we then obtain 
[B j ,Bi]{t) = [V , L J ]{t(t)) = f* (t) H{s)ds = I For 
i ^ j, we have , B l ]{i) = f* (t) H(s) d[B j , B%s) = 0. 
Thus, Levy's Theorem [55 implies that B(t) = {B j (t)) 
is a d-dimensional standard Wiener process. 

Finally, using the definitions of Y, E, and F J , we 
find [H 



Y(t) 



Ki) 



E(s) ds 



t(t) 



F J (s)dB j (s) 



E(t(s)) 



ds 



o H(t(s)) J y/H(t(S)) 
E(s)ds + [ Fji^dB^s), 



dB^i 



(13) 



which is just the SDE (|12c|) written in integral notation. 

Conclusions.- The above discussion shows how rela- 
tivistic Langevin equations can be Lorentz transformed 
and reparameterized within a common framework. Thus, 
mathematically, th e sp ecial relativistic Langevin the- 



ory ll3|) UA, UM, lis LLZ( is now as complete as the clas- 



sical theories of nonrelativistic Brownian motions and 
deterministic relativistic motions, respectively, both of 
which are included as special limit cases. From a physics 
point of view, the most remarkable observation con- 
sists in the fact that the r-parameterized Brownian mo- 
tion converges to a modified Juttner function fiol ] if 
the corresponding ^-parameterized process converges to 
a Juttner function [381 ]. This illustrates that it is nec- 
essary to distinguish different notions of 'stationarity' in 
special relativity. While the i-parameterization appears 
more natural when describing diffusion pro cesses from 

the viewpoint of an external observer m 111 lis mm 



23, 24 



25j, the r-parameterization is more convenient 
when extending the above theory to include particle cre- 
ation/annihilation processes, because a particle's lifetime 
is typically quantified in terms of its proper-time r. Last 
but not least, the proper-time parameterization paves the 
way toward generalizing the above concepts to general 
relativity. 
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